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Abstract. The spin modulated gravitational wave signals, which we shall call 
smirches, emitted by stellar mass black holes tumbling and inspiralling into massive 
black holes have extremely complicated shapes. Tracking these signals with the aid of 
pattern matching techniques, such as Wiener filtering, is likely to be computationally 
an impossible exercise. In this article we propose using a mixture of optimal and 
non-optimal methods to create a search hierarchy to ease the computational burden. 
Furthermore, by employing the method of principal components (also known as singular 
value decomposition) we explicitly demonstrate that the effective dimensionality of the 
search parameter space of smirches is likely to be just three or four, much smaller than 
what has hitherto been thought to be about nine or ten. This result, based on a limited 
study of the parameter space, should be confirmed by a more exhaustive study over 
the parameter space as well as Monte-Carlo simulations to test the predictions made 
in this paper. 
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1. Introduction 

A stellar mass compact object, either a neutron star or a black hole, in tight orbit around 
a massive Kerr black hole will tumble around the massive companion due to strong 
spin-orbit coupling. For small mass ratios the stellar mass companion slowly inspirals 
into the massive hole, the duration of inspiral increasing inversely as the mass ratio. 
Gravitational waves emitted in the process of adiabatic inspiral encode the structure 
of the spacetime geometry around the hole thereby allowing us to test the uniqueness 
theorems of black holes spacetimes [1]. However, decoding the information carried by 
the waves would require an accurate knowledge of how the emitted waves depend on 
the parameters of the system and fitting the data with the expected signal. There are 
problems with both of these aspects: On the one hand we do not currently have an 
accurate knowledge of the signal although over the next decade the situation is likely 
to have improved. On the other hand, the problem of fitting is unlikely to go away 
as the number template waveforms required to decode the useful information might 
be formidably large. This is because the dimensionality of the parameter space is quite 
large. Indeed, in the generic case of a binary consisting of two spinning black holes there 
are potentially seventeen parameters to fit: The direction and distance to the source, 
two masses, three vectors corresponding to the initial values of the orbital angular 
momentum and two spin angular momenta, the eccentricity of the orbit, the instant 
when the orbital frequency reaches a fiducial value and the phase of the signal at that 
instant. 

Even though the time-structure of the signal depends on seventeen parameters, 
varying different parameters may not always produce distinct waveforms. In other 
words, the effective dimensionality of the parameter space for the purpose of fitting 
might be a lot smaller. For instance, we know in the case of non-spinning binaries that 
out of potentially eight parameters - namely, the direction to the source and its distance 
from the Earth, two masses, an initial magnitude of the angular momentum, a fiducial 
time and phase when the signal reaches a certain frequency - there is effectively only 
one parameter with which to fit the data. This is the so-called chirp mass, which is a 
certain combination of the two masses of the component bodies. The detector response 
has no time-varying dependence on the source direction since the signal lasts for a 
very short duration (minutes, at most) and therefore doesn't suffer any modulation 
due to the motion of the antenna that could have been useful in deciphering the 
source direction; the initial magnitude of the angular momentum only sets an initial 
frequency of observation that is chosen by hand; the fiducial time and phase can be 
easily determined by the use of the fast Fourier transform and analytically, respectively; 
and the phase of the signal, which is what is important in accomplishing a good fit, 
depends, at the dominant post-Newtonian order, only on the chirp mass. Although the 
degeneracy in component masses is resolved when higher order post-Newtonian terms 
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are included in the phasing, it turns out the number of 'independent' shapes doesn't 
increase significantly while including a second mass parameter. 

An analysis similar to the non-spinning case is required in the case of spinning 
black hole binaries. In this paper we discuss the method of principal components [2] 
that can be used to determine the effective dimensionality of the search space. We will 
then apply it to the case of a test mass falling into a massive black hole. We sample 
the parameter space at a few crucial points to get an idea of the number of principal 
components on which the signal shape depends, thereby obtaining an estimate of how 
large the parameter space could be. We begin with a rough counting of the number of 
template waveforms that may be required in a fitting problem. This is really an upper 
limit on the number of waveforms, the actual number might be a lot less. 

Since we expect the number of templates to be too large it will not be possible to 
employ matched filtering on an entire year's worth of data. Rather, the idea is to employ 
a technique, such as time-frequency analysis, that makes little use of the signal shape in 
capturing a signal, albeit with a lower detection probability. Such an analysis might help 
to first get a rough idea about the signal parameters. One could then employ matched 
filtering on, say a week's worth of data and refine the estimation of parameters. Based 
on these refined parameters one finally employs the search on a year's worth of data but 
in a small region of the parameter space. In other words, we are advocating a multi-step 
hierarchical search involving different analyses techniques. In such a scheme it would 
only be necessary to estimate the number of templates required in the most compute 
intensive stage of the search hierarchy, namely in the first stage of matched filtering. 
It is best to begin the search from a point where the parameters can be nailed down 
most easily, namely by employing matched filtering on a week's worth of data. (At this 
stage a signal to noise ratio of about 5-10 can be expected for mergers taking place at a 
distance of about a few Gpc). This assumption somewhat simplifies our calculation as 
we can neglect the motion of LISA relative to the source during this period and assume 
that the detector is fixed relative to the source. Strictly speaking, this is not a valid 
assumption but it is unlikely to change the qualitative conclusions reported in this work. 

2. Counting the number of templates — A rough estimate 

To count the number of templates we must first understand how the signal evolves as 
a function of time. To this end we shall use the best post-Newtonian order currently 
known but ignore eccentricity. It turns out that at this order the signal already has many 
modulations including spin-orbit and spin-spin couplings. Therefore, the waveform 
should be quite suitable for counting template numbers although it may not be good 
enough for the purpose of fitting. 
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2.1. Evolution of the orbit 

Following Ref. [3] we shall use a coordinate system adapted to LISA. In the post- 
Newtonian approximation the evolution of a binary system comprising of two bodies 
of masses mi and m 2 , spins Si = (s 1 , 9 1 , ipi) and S 2 = (s 2 , 9 2 , ^2), orbital angular 
momentum L = (L, 9q, ipo), is governed by a set of differential equations given by: 



Ami + 3m 2 3S 2 -L\g /4m 2 + 3mi 3 Si • L N 
2mim 3 2 L 2 m 3 J 1 y 2m 2 m 3 2 L 2 m 3 / 



3 V^ T 7 m 



Si = 



'4mi + 3m 2 3 S 2 ■ L\ ^ ^ ^ S 2 x Si 

2mim 3 2 L 2 m 3 / 1 2m 3 

' 4m 2 + 3mi _ 3 Si ■ L \ L x g + g_i x S 2 

2m 2 m 3 2 L 2 m 3 / 2 2m 3 



(2) 



(3) 



Here v = (irmf) 1 / 3 , where / is the gravitational wave frequency, is the post-Newtonian 
velocity parameter, an overdot denotes the time- derivative, m = mi + m 2 is the total 
mass, and r\ = mim 2 /m 2 is the (symmetric) mass ratio. In the evolution of the orbital 
angular momentum we have included the lowest order dissipative term [ the second term 
containing v 7 in Equation (1) ] while keeping the non-dissipative modulation effects 
caused by spin-orbit and spin-spin couplings [the first term within square brackets 
containing v 6 in Equation (1) ]. The spins evolve non-dissipatively but their orientations 
change due to spin-orbit and spin-spin couplings. Though the non-dissipative terms are 
not responsible for gravitational wave emission, and therefore do not shrink the orbit, 
they cause to precess the orbit. 

In the absence of spins the antenna observes the same polarization at all times; 
the amplitude and frequency of the signal both increase monotonically, giving rise to a 
chirping signal. Precession of the orbit and spins cause modulations in the amplitude 
and phase of the signal and smear the signal's energy spectrum over a wide band. The 
smearing is particularly dramatic when one of the bodies is much lighter than the other 
such as the inspiral of a stellar mass compact object into a massive black hole. In this 
case the orbit would undergo many precessions before the small object plunges into the 
hole smearing the signal into hundreds of frequency bins as opposed to a few frequency 
bins if the signal were a simple chirp. For this reason we shall call the spin modulated 
(sm) chirp, a smirch (an anagram of sm and chir). 

2.2. The waveform 

The strain h(t) produced by a smirch in LISA is given by 

h(t) = -A(t)cos[$(t)+<p(t)], (4) 

where A(t) is the precession- modulated amplitude of the signal, $(£) is its post- 
Newtonian carrier phase that increases monotonically and ip(t) is the polarization phase 
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caused by the changing polarization of the wave relative to the antenna. (We have 
neglected the Thomas precession of the orbit which induces additional, but small, 
corrections in the phase.) For a source with position vector N = (D, 9s, <ps) the 
amplitude is given by, 



A(t) 



2r\mv 1 
D 



(l + (L • N) 2 ) 2 F+(9s, Vs , + 4 (L • N) 2 F 2 x (9 s , <p s , i/>) 



1/2 



(5) 



Here L = L/L, N = N/D (D is the distance to the source) and ip(t) is the precession- 
modulated polarization angle given by 

L(t) - Z- (L(t) -N)(z-N) 



tan0>(t) 



(6) 



N • (L(t) x z) 

Also, F + and F x are the antenna beam pattern functions are given by 

F + (9 S , ips, VO — ^ (l + cos2 #s) cos 205 cos 2-0 — cos 9 S sin 20 5 sin 2-0, (7) 

F x (9s, 0s, V') — \ (l + cos2 cos 205 sin 2-0 + cos #s sin 20 5 sin 2-0. (8) 



Next, the polarization phase ip(t) is 

2L(t)-NF x (9 s ,<Ps,'4>) 



tan 



(9) 



l + (L(t)-N) 2 jF + (^,^,0;) 

And finally, for the carrier phase we use the post- Newtonian expression, but without the 
spin-orbit and spin-spin couplings. These spin couplings modify the carrier phase by 
amounts much smaller than the post-Newtonian effects so that neglecting them makes 
no appreciable difference to our main conclusions. To second post-Newtonian order the 
carrier phase is given by 



*(*) = 



-2 

^05 



1 + 



/3715 55 



8064 96 



3tt 
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V 14450688 258048 



1855^ 
2048 5 



v 2 )e 4 



.(10) 



where 6 = [r)(tc — t) / (5m)] 1 / 8 , to being the time at which the two stars merge together 
and the gravitational wave frequency formally diverges. 



2.3. A rough count of templates 

From the foregoing equations we can get an idea of the number of parameters that may 
be needed in a search. The system we are interested in is that of a stellar mass black hole 
falling into a supermassive black hole wherein the mass ratio could be as low as 10~ 5 . 
Consequently, the spin of the stellar mass hole, being proportional to the square of its 
mass, would be negligible compared to the spin of the massive hole. Further, we don't 
have to specifically search for the distance to the binary nor for the fiducial time and 
phase when the frequency of the emitted gravitational wave reaches a certain value as 
they can all be measured without incurring significant computational costs. Moreover, 
since we are focussing on quasi-circular orbits (i.e. zero eccentricity) and only those 
systems that merge within a week of observation (i.e. after about 1000 orbital time 



Templates for Smirches 



6 



scales) the number of parameters in our search is down from 17 to 9. These are the 
initial direction of the orbital angular momentum (2 angles), the spin of the massive 
black hole (1 magnitude and 2 angles), the masses of the two bodies (2 numbers) and 
the direction to the source (2 angles). 

As mentioned in the introduction we don't expect both masses to be important but 
only a combination - the chirp mass M. = rf^m. The azimuth angles are normally not 
important in a search but co-latitudes would lead to modulations that cannot be easily 
searched for. Consequently, out of the six angles we expect only three to be important. 
(Note that due to different orientations of the binary and the antenna these angles get 
mixed up and therefore azimuth and co-latitudes are not well-defined. However, we 
expect certain combinations of these angles to play the roles of azimuth and co-latitude. 
Our claim is that three of those combinations would be important in a search and the 
other three not so important.) In addition, we expect the spin magnitude of the hole to 
be a search parameter although this would be important only when s± is close to 1. As a 
result the number of search parameters would be five or less. Precisely how many would 
have to be determined by a more rigorous and quantitative analysis. Assuming that 
we have to search for five parameters in all and that each parameter requires as many 
templates as the number of cycles we get a rough count of (10 3 ) 5 = 10 15 templates. We 
shall see below that this is a rather large over estimate. The number of independent 
parameters is possibly about three leading to a significant reduction in the number of 
templates. 



3. Covariance matrix, correlation coefficients and principal components 

Any measurement process is disturbed by a background perturbation or noise which 
causes a systematic and/or random errors in the measured quantities. In general, 
the errors get smaller at greater signal-to-noise ratios. We review the theory behind 
estimating errors involved in a measurement process. 

For a signal buried in a Gaussian background and detected using matched filtering 
the covariance matrix Hki - a symmetric matrix whose diagonal elements are the 
variances in the measurement of the various parameters and the off-diagonal elements 
are the covariances between them. The covariance matrix is itself the inverse of the 
information matrix, or the metric, g mn . Let us suppose the signal h(t; A&) depends on p 
parameters A*,, k — 1, . . . ,p. The information matrix is defined by 

g m n = (d m h , d n h) , S mra = g' 1 , (11) 

L J mn 

where for any two functions a(t) and b(t) (and Fourier transforms a(f) and b(f)) their 
scalar product (a , b) is defined as: 

(a,b) =2 / /high ^T [«(/)&*(/) + r(/)&(/)] , (12) 

-'/low Jh{j) L J 

where Sh(f) is the (real) one-sided noise spectral density of the antenna and /i ow and 
/high are some appropriately defined lower and upper frequency cutoffs depending on 
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the signal and/or the antenna response. 

Starting from the covariance matrix one defines the matrix of correlation coefficients 
Cu in the following manner: 



As is obvious the diagonal element Ckk of the matrix of correlation coefficients is the 
(one-sigma) standard deviation in the measurement of the parameter A&, and the kl-oS- 
diagonal element Cm is the normalized covariance between the parameters A& and A/. 
The correlation coefficients C^i take values in the range [—1, 1]. A correlation coefficient 
Cki close to 1 (or —1) indicates that the parameters A& and A/ are perfectly correlated 
(or anti-correlated) with each other. If on the other hand the value is close to zero then 
there is little correlation between the two parameters. When two parameters have a 
correlation coefficient close to ±1 then varying one of the parameters is equivalent to 
varying the other. In other words one doesn't produce two independent signal shapes by 
varying both the parameters, varying only one of them should be sufficient. However, 
when two parameters are not correlated then varying one of them produces a distinct 
shape as compared to varying the other. Therefore, we can look at the matrix of 
correlation coefficients to determine how many parameters are really independent. 

Finally, we can diagonalize the covariance matrix to find the principal components. 
On diagonalization one has (linearly) transformed from the parameter set A& to a new set 
Hk effectively removing in the process all the covariances between the parameters. The 
principal components of the signal are those transformed parameters that correspond 
to the largest diagonal elements. Since we have (by construction) a symmetric matrix 
the eigenvalues are all real. Moreover, since the covariance matrix is the inverse of the 
metric the principal components are the inverses of the eigenvalues of the metric. Before 
we embark upon a discussion of the results let us look at two simple examples. 

Our first example is a simple mathematical example: Consider a two-dimensional 
matrix C k i. Suppose Cu = C 22 = 1 and C 12 = C 21 = e < 1. The eigenvalues A± of such 
a matrix are A± = 1 ± e. If e is close to zero then both the eigenvalues are roughly the 
same and both the parameters are equally important. If, however, |e| is close to 1 then 
one of the eigenvalues is nearly zero making the corresponding transformed parameter 
unimportant. Thus, closer the correlation coefficient is to ±1, smaller is the significance 
of the parameter. 

As a second example let us consider a sinusoidal signal h(t; A, to, (po) = A cos[u(t — 
t ) + ipo\. Here the parameters A and t are uncorrelated and so are A and f , while t 
and ipo are perfectly anti-correlated. Increasing ip is equivalent to decreasing t but the 
signal has really only two independent parameters leading to two principal components. 





if k^l. 



(13) 
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Table 1. Correlation coefficients and errors in the estimation of parameters for a 
smirch corresponding to the following parameters: m = 1O 6 M0, r\ = 10~ 3 , a\ — 
0.9, 6i = 7r/4, tpi = it/2, 9l = 7r/5, y>L = tt/3, 6*s = 7r/8, ips = n and n cyc = 910 





m 


V 


ai 


0i 




0l 


PL 


0s 


<Ps 


to 


0.0033 


















7/ 


0.9592 


0.0023 
















di 


0.9336 


0.8775 


0.0035 














^1 


0.1516 


-0.1776 


-0.1301 


0.3270 












¥>1 


-0.0275 


0.0326 


0.0262 


0.7073 


1.8633 










#L 


0.0771 


-0.0840 


-0.0973 


0.7127 


0.3228 


0.5821 










-0.0449 


0.0526 


0.0472 


0.6975 


0.9984 


0.2855 


2.036 






0s 


-0.0002 


-0.00341 


-0.0101 


-0.2898 


0.0269 


0.8667 


0.0690 


0.8522 




<PS 


-0.0150 


0.0194 


0.0197 


0.8135 


0.9055 


0.6809 


0.8875 


-0.3925 


1.3879 



4. Correlation coefficients for smirches 

We now apply the formalism developed in the previous section to the gravitational 
wave signals produced by stellar mass black holes falling into supermassive black holes. 
We first compute the covariance matrix and then diagonalize it to obtain the principal 
components. 

Let us first examine the matrix of correlation coefficients. In Tables 1-3 we have 
displayed the correlation coefficients of the smirches for three different values of the 
mass ratio, rj = 10 -3 , 10 -4 , 10~ 5 . We have also changed the orientations of the spin, 
orbital angular momentum and source direction from one example to the other to get 
an idea of how the correlation coefficients depend on the angular parameters. Tables 1-3 
correspond to cases A, B and C, respectively; in all cases the total mass m = 10 6 M Q , 
ai = 0.9 and the initial orientation of the orbit is the same (Olj^l) = (tt/2, 7r/3). 

r] 6>i ip x 9 S ip s 

Case A: 1(T 3 vr/4 tt/2 tt/8 tt 

Case B: 1(T 4 tt/3 tt/4 tt/2 tt/3 

Case C: 1(T 5 tt/4 tt/2 tt/2 tt/3 

The results reported here should be taken as a hint on what the general results are 
going to be but we cannot be absolutely certain about any of the results. We plan to 
report a fuller treatment of the problem with exhaustive Monte-Carlo simulations to 
validate our results. 

The diagonal elements in these Tables are the expected standard deviations in 
the measurement of the corresponding parameters when the (amplitude) signal-to-noise 
ratio is 10. In the case of m, i] and s\ we have quoted the relative error (namely, Uk/Xk 
where a k is the expected standard deviation and A fc is the true value of the parameter) 
while the error in angles is the absolute error expected in their measurements. Since the 
signal-to-noise ratio is proportional to y/r]n cyc we have chosen the quantity y/i]n cyc to 
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Table 2. Correlation coefficients and errors in the estimation of parameters for a 
smirch corresponding to the following parameters: m = 10 6 M Q , rj = 10~ 4 , a\ — 







0.9, 01 = 7T 


A <pi = 


3tt/4, e L = 


n/5,tp L = 


tt/3, s = 


tt/2, fs = 


7r/3 and 


n cyc — 






9.1 x 10 3 
















k 


m 


n 




0i 


t/?1 


L 


T 

r Aj 


0s 


(J) a 


m 


0.0003 


















V 


0.9449 


0.0002 
















ai 


0.9782 


0.9358 


0.0003 














Oi 


0.0057 


-0.0068 


-0.0184 


0.5922 












<fii 


0.0094 


-0.0113 


-0.0091 


0.6862 


0.8631 










Ol 


0.0087 


-0.0160 


-0.0169 


0.9847 


0.7009 


1.0201 








PL 


0.0200 


-0.0263 


-0.0331 


0.8773 


0.9425 


0.8675 


1.2781 






Os 


0.0096 


-0.0140 


0.0014 


0.9954 


0.6963 


0.9935 


0.8761 


1.0203 




<PS 


-0.0082 


0.0100 


0.0106 


-0.1538 


0.5862 


0.1540 


0.3332 


0.1550 


0.6145 



be roughly the same in the three examples^. For this reason, though our sources are all 
at the same distance they produce similar signal-to-noise ratios. In reality they have to 
be integrated over different time scales, the duration being proportional to 77 -1 . Finally, 
note that though we have quoted a total mass (essentially the mass of the supermassive 
black hole) of 1O 6 M the results of this analysis hold good for other masses as well 
provided the mass ratios are the same. The spectral density of noise really doesn't vary 
too much over the duration of a signal so that the results discussed here are also valid 
for ground-based detectors with some minor modifications. 

The values of the correlation coefficients whose absolute value is greater than 0.5 is 
shown in bold-face to easily identify which parameters are correlated and which aren't. A 
correlation coefficient of 0.5 does not mean that one of the parameters can be completely 
ignored. However, it does mean that the two parameters are not completely independent 
and we can save some computational cost by varying a suitable combination of the two. 
If the correlation is close to 1, however, we could save a lot in computational costs. 

In respect of correlation coefficients the parameters fall into two categories: the 
angular parameters describing the orientation of the black hole's spin, orbital angular 
momentum and source location fall into one category and the non-angular parameters 
corresponding to the magnitude of the massive black hole's spin and the two masses fall 
into another category. The parameters within each set are correlated with one another 

% Indeed, as a rule of thumb all smirches with the same total mass m but different 7/(<C 1), that begin 
at the same frequency and coalesce within the LISA band will have the same signal-to-noise ratio. This 
is because, if the starting frequency is the same, then as r\ is lowered by a certain factor the number of 
cycles, and the duration of the signal, goes up by the same factor so that r\n cyc is the same. Even though 
the duration of the signal changes quite a lot the frequency content of the signal doesn't change so much 
during this interval as a result of which the signal-to-noise ratio is roughly invariant. For ground-based 
sources such a rule of thumb doesn't apply since these sources span the entire detector band and the 
signal-to-noise ratio depends not only on the factor r)n cyc but also on where in the frequency band 
those cycles are located. 
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Table 3. Correlation coefficients and errors in the estimation of parameters for a 
smirch corresponding to the following parameters: m = 10 6 M Q , r\ = 10~ 5 , a\ — 







0.9, 0i = 7i 


'/4, <pi = 


tt/2, l = 


n/5,(pL = 


tt/3, s = 


t/2, ip s = 


7r/3 and 


^-cyc — 






9.1 x 10 4 
















k 


m 


V 




0i 




0l 




0s 


(fig 


m 


2.5e - 5 


















V 


0.9370 


1.4e-5 
















ai 


0.8804 


0.7873 


2.1 e- 5 














0i 


-0.0315 


0.0435 


0.0370 


0.4220 












<Pi 


-0.0003 


0.0018 


-0.0182 


0.8500 


1.1942 










Ol 


-0.0124 


0.0188 


-0.0057 


0.9781 


0.8812 


0.9363 








PL 


-0.0150 


0.0219 


-0.0047 


0.7517 


0.9723 


0.8242 


1.1773 






0s 


0.0330 


-0.0490 


-0.0097 


0.9186 


0.8790 


0.9733 


0.8637 


0.9398 






0.0040 


-0.0063 


-0.0010 


0.1735 


0.2853 


0.1819 


0.2969 


0.1847 


0.5546 



but there is little correlation between parameters from the two categories. 

Let us return to the examples given in Tables 1-3. The parameters m, rj and si 
are so strongly correlated that only one of these three parameters suffices in a search. 
Similarly looking at the correlation between the angular parameters we see that just two 
or three of the angular parameters might be good enough in a search. In the case of the 
masses and spin magnitude the correlation coefficients remain roughly the same in the 
three different examples. However, this is not so in the case of angular parameters. For 
the sake of defmiteness we declare a parameter to be independent, and would therefore 
be needed in a search, if its correlation coefficient with the other parameters is smaller 
than 0.85. Comparing Tables 1-3 we find that in the first example (Table 1) 9 1 , ipi, and 
L are the three independent angular parameters while in the second example (Table 2) 
6i, <fi, and Lps are the three independent parameters and in the third example 6\ and 
ips are the independent parameters. If the independent parameters had remained the 
same then the search algorithm would be a lot simpler; we could determine forever the 
independent parameters and carry out a search in that space irrespective of where in 
the parameter space we are. However, now the independent parameters for our search 
(in a small region of the parameter space) depends on where in the parameter space we 
are. One should measure the matrix of correlation coefficients at a given point in the 
parameter space before embarking upon laying a grid of templates. 

5. Principal components for smirches 

Next we turn to the principal components of the spin modulated chirp. Recollect that the 
principal components are nothing but the eigenvalues of the covariance matrix. In Table 
4 we have shown the principal components for the cases discussed in Tables 1-3. We have 
displayed only those components that are a significant fraction of the largest component. 
(A rigorous treatment of the problem should deal with dimensionless parameters so that 
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Table 4. Principal components for the three cases discussed in Tables 1-3. Only 
components greater than 1CP 2 (chosen in an ad hoc manner) are listed. Cases A-C 
correspond to the systems discussed in Tables 1-3 (see text) . 



Case A Case B Case C 



13 


445 


398 


9.3 


4.3 


4.3 


0.19 


0.025 


0.11 


0.15 




0.016 


0.12 






0.050 







in the process of diagonalization the transformation of the parameters is meaningful. 
We defer such a treatment to a later paper.) We note that the number of principal 
components is much less than the number of parameters although the precise number 
of principal components depends on the point in the parameter space. Our preliminary 
conclusion is that the problem of catching smirches boils down to searching for three, 
or at most four, parameters. This means that we might need a lot fewer templates than 
was originally conceived. After all the problem may really not be as serious as we had 
originally thought. 

6. Conclusions 

Gravitational waves emitted by stellar mass black holes falling into a supermassive black 
hole encode a lot of information about the structure of the spacetime geometry around 
the massive Kerr black hole. Decoding this information requires not only accurate 
modelling of the emitted waveform but algorithms that can efficiently dig out and 
discriminate these signals. Since the parameter space corresponding to such systems 
is quite large it has been thought that matched filtering, which is not only effective in 
detecting signals but also deciphering their information content, would be too expensive 
in a search for these signals. In this paper we have explored the signal manifold 
corresponding to spin modulated chirps, or smirches as we would like to call them, 
using the covariance matrix and the method of principal components. We conclude 
that the problem might be tractable after all since the effective dimensionality of the 
parameter space might be a lot smaller than what has been thought before. 

The signal emitted by a small black hole in quasi-circular orbit around a massive 
Kerr black hole is, in general, characterized by nine search parameters - the two masses, 
the magnitude of the black hole's spin and its orientation, initial orientation of the orbit 
and the direction to the source. Our analysis has shown strong correlations between the 
two masses and the spin magnitude as well as between the different angular parameters. 
We have found strong indications that matched filtering a short duration smirch (1,000 
to 100,000 cycles) might require a search in only only 3-, or at most 4-dimensional 
parameter space. This means that one might be able to use matched filtering to fit the 
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strongest smirches in LISA's data and dig deeper for weaker signals. Before we can be 
completely optimistic about the result, a Monte-Carlo simulation should be performed 
to see if the overlap of arbitrary smirches can be larger than a certain minimal match 
(say 50%) when maximized over just three parameters. 
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